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ABSTRACT 



We have produced deep radio maps of the double quasar 0957+561 from 
multiple-epoch VLA observations. To achieve high sensitivity to extended 
structure we have re- reduced the best available 1.6 GHz observations and 
have combined 5 GHz data from multiple array configurations. Regions of 
faint emission approximately 15" north and south of the radio source G are 
probably lobes associated with the lensing galaxy. An arc 5" to the east of G 
may be a stretched image of emission in the background quasar's environment. 
1.4" southwest of G we detect a source that we interpret as an image of emission 
from the quasar's western lobe, which could provide a constraint on the slope 
of the gravitational potential in the central region of the lens. We explore the 
consequences of these new constraints with simple lens models of the system. 



Subject headings: gravitational lensing — quasars: individual (0957+561) 



1. Introduction 

Astrophysicists have anticipated the use of gravitational lensing as an observational 
tool for 60 years ( [Zwicky 1937| ; |Schneider, Ehlers, fc Falco 1992j ), and in the case of the 
double quasar 0957+561 ( |Walsh, Carswell, fc Weymann 1979] ), after 18 years of study 
the promise is closest to fulfillment. If one knew the details of the gravitational- lensing 
potential, and the time delay among the images of flux-variable components, one could 
make an estimate, albeit cosmology-dependent, of Hubble's constant H ([Refsdal 1964j) . 

Efforts to measure the time delay in this system have converged recently (417+3, 
Kundic et al. 1997] 420+13, [Haarsma 1997|) . However, models of the lensing potential have 
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been less well constrained ([Falco et al. 199 1 ; Kochanek 1991 ; Grogin fc Narayan 1996a 



1996b) despite detailed observations of the cluster of galaxies providing the lensing mass 
( |Young et al. 198T| ; |Angonin-Willaime, Soucail, fc Vanderriest 1994] ; [b'ischer et al. 1997| ). 
In an effort to produce a definitive radio map of the object we undertook to re-reduce 
VLAQ data gathered by the M.I.T. group, discovering several new features in the field 
( |Avruch et al. 1993|) . In this letter we present improved maps and identify features that 
may be useful as model constraints. 



2. Observations 

The data sets from which results presented in this letter were computed are listed 
in Table 1. We first mapped a low resolution A6cm data set to identify any sources 
in the primary beam whose side lobes would contaminate the field of interest; during 
self-calibration of other data sets this emission was taken into account. For sensitivity 
to low surface brightness features we chose the best extant A18cm observation. We have 
also combined five A6cm data sets from array configurations A, B, and C to achieve more 
complete (u, v) coverage and compensate for the reduced brightness at A6cm compared 
to A18cm. The fluxes of A and B were roughly constant at these epochs. Using the 
MPS software, the data sets were independently mapped and self-calibrated following 
standard VLA reduction procedures (|Fomalont fc Perley 1989|) . Each data set was phase 



self-calibrated several times, followed by a single amplitude self-calibration, provided that 
it reduced the map noise. The individual data sets were then co-added in AIPS and the 
combined data were mapped and self-calibrated as above. To produce source-subtracted 
images, we use the model for compact emission that the deconvolution algorithm creates in 
the form of CLEAN components, subtracting the model source from the visibility plane and 
remapping. 

To the north and south of 0957+561 we have detected lobes of emission, separated by 
about 30". The northern lobe (N) is more compact, with a A6cm flux of about 840 /iJy and 
spectral index a 1 ^™ ~ —1.0 (S oc v a ). The southern lobe (S) is extended, with a total 
flux of about 1000 /iJy, a 1 ^™ ~ —0.7. These lobes may be associated with radio galaxy G, 
or with the lensed quasar in the background. To the east of the quasar images A and B 
we have detected an arc of emission (Rl). The arc is clearly resolved tangentially, with a 
peak flux of 1.27 mJy beam -1 at A18cm and spectral index a 1 ^™ ~ —0.8. In Figure 1 we 



2 The VLA is part of the National Radio Astronomy Observatory, which is operated by Associated 
Universities, Inc. under co-operative agreement with the National Science Foundation. 



- 3- 



present a radio map of these features; B has been subtracted from the image in the manner 
described above. 

Galaxy G, the dominant contributor to the lensing potential, is definitely extended to 
the east, southwest, and northwest. To better view the structure near G, we subtracted 
from the multi-epoch (u, v) data all emission associated with the B quasar image, the 
BN component ([Roberts et al. 1985| ), and G. These structures were identified by directly 



inspecting the CLEAN components from the multi-epoch map. In Figure 2 the extension of 
G to the east we name GE, to the northwest GN, to the northeast GNE, and the brightest 
component of the arc-like structure to the southwest of G we call R2. Table 2 presents the 
positions and fluxes for these new components. 

We are confident that these features are real. N, S, and Rl have been confirmed 
with detections by Harvanek et al. (1996); Rl and perhaps GN have also been confirmed 
by Porcas et al. (1996). The fainter features GE, GN, GNE, and R2 are visible in every 
individual, reduced data set with sufficient resolution and sensitivity, so it is unlikely that 
they are artifacts of calibration or deconvolution. On the other hand, detailed substructure 
such as the double peaks of GNE is not significant, because with extended sources CLEAN 
produces spurious peaks on that scale ( |Briggs 1995| ). 



3. Discussion 

To illustrate our interpretation of these new VLA components, we used the LENSMOD 
software QLehar et al. 1993] ) to model the lensing mass with a softened power-law potential 
( [Blandford fc Kochanek 1987] ). The model parameters were: the lens position (Act, 
AS), the critical radius (b), a core radius (9 C ), the power index P (P = 1 is isothermal, 
while P = 2 is a Hubble profile), the isodensity ellipticity (e = 1 - mmorax ! s ), and the 

r 11 j sr j \ major axis ' ' 

major axis orientation (</>). As constraints we used the new HST quasar and Gl positions 
( Bernstein et al. 1997| ) and required that the quasar images have a magnification ratio of 



0.75±0.02 ( [Schild fc Smith 1990| ). We required that any third image of the quasar near G 
be at least 30 times fainter than B (as a la limit). We also added constraints from the new 
HST "blobs" and "arc." We required that blob2 and blob3 be images of each other, and 
that the two knots in the arc share a common source. Note that the HST arc is probably 
caused by the eastern end of the same object that gives rise to blobs 2 and 3, and this 
could be used to further constrain lens models. To account for the possibility that the 
HST objects are at a different redshift than the quasar, we added a uniform scale factor 
Q2 to the deflection angles for those components, as an extra model parameter. The lens 
model parameters were varied until the source plane position and magnitude differences 
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for each pair were minimized, with a resultant reduced x 2 f° r the fit of 1.1. The best fit 
model parameters are given in Table 3, with uncertainties determined by varying the model 
parameters until the reduced x 2 increased by 1. Note that the Q 2 range corresponds to HST 
component redshifts of zhst ~ 1-3 ±0.1 for an Q = 1 cosmology, which is fully consistent 
with the quasar and HST objects being at the same redshift. Figure 3 shows the best fit 
model for Q 2 — 1, with components added to show the modeled radio emission. We do not 
attempt to account for the VLBI structures QGarrett et al. 1994] ) in this model, and thus 



make no claims about the time delay or Hubble's constant based upon our model. 

We interpret the component GE as the counter-image to the low surface brightness 
tail of the quasar's western radio lobe E. GE's peak surface brightness and spectral index 
( al 6cm ~ —1.0) matches that of component E's northeastern extension, so the brighter parts 
of the lobe are not multiply imaged. The Bernstein et al. (1996) HST blobs 2 and 3, almost 
certainly multiple images of a background object, are very close to the positions of GE 
and the northeast end of E; therefore we expect an image of E near where we have found 
GE. Because not all of E is multiply imaged, the detailed structure of GE can yield strong 
constraints on the central region of the lens: either the mass distribution is non-singular, in 
which case GE comprises two merging images of the eastern end of E, or, if the mass has 
a central singularity, GE will have a sharp cusp at its western end. High resolution radio 
observations of GE may be able to distinguish these two possibilities, or at least determine 
an upper limit on the size of the central mass concentration in G. This is also important 
because, for a given lens mass, the potential near the quasar B image is generally deeper for 
singular models, yielding a longer predicted time delay and thus a lower Ho estimate. 

The arc-like feature Rl may be a stretched image of background emission. As there 
is no clear counterpart to the west of G, it is unlikely to be multiply imaged. If the 
background source is circular, the axial ratio of Rl yields a lower limit of about 5 for its 
magnification. Jones et al. (1993), in Einstein HRI data, have detected an apparent x-ray 
arc about 3" northwest of Rl. The positions are formally consistent, but seem unlikely to 
be coincident judging from the relative positions of A and B. An association is not ruled 
out, however. The authors claim the extended x-ray source is an image of thermal emission 
from a cooling flow in the cluster hosting the lensed quasar at z = 1.41. There are examples 
of diffuse non-thermal radio emission associated with x-ray-emitting clusters ( Peiss et al 



1997|) , and in this case the lensing magnification may have helped to make it observable. 
Of course this emission could be foreground; if G has radio lobes, N and S, it could as well 
have jets. Rl might be back flow from the lobe S, and GN might be a faint jet feeding the 
lobe N. GE is well explained as an image of the quasar's E lobe, but it's not impossible for 
it to be the counter-jet of GN, feeding lobe S. 
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The features R2, GN, and GNE are not readily explained by a lensing hypothesis. R2 
is in the position of the western half of the HST arc, but all the models we have investigated 
would produce an eastward extension of this arc which is not detected. We could appeal, 
ad hoc, to source size and spectral index morphology causing the image to be unobservable. 
The component GN should have a brighter image 5" south of G, which is not seen, though 
we could make the same appeal and note that it might be difficult to separate visually from 
S. GNE should also have a counter-image to the south of G, which is not seen. However, 
given the interpretation of Rl as lensed, and the faintness of these features, it is not ruled 
out that at least some of the emission is due to structure in the background quasar's 
environment. 

N and S are certainly not multiply imaged, but whether they are foreground or 
background is less clear. They could be the radio lobes of the galaxy G. At the lens redshift 
(z = 0.36, and assuming Q — 1, h — 0.75) N and S would have a (projected) proper 
separation of 120 kpc, and luminosities at 178 MHz of about 10 24 WHz -1 , typical values for 
low power, limb darkened radio galaxies. The optical classification of G as a cD galaxy, and 
the fact that N and S are aligned within 30° of its optical minor axis are also consistent 
( Milcy 1980| ). N and S might be old lobes of the background quasar, in which case the 
numbers are 170 kpc and 10 26 WHz _1 , more appropriate for powerful, limb brightened 
sources. If N and S are associated with the quasar, the relatively small lobe separation (56 
kpc) and the high core-to-lobe flux ratio (R = 0.22) suggest that the jet axis is moderately 
inclined towards the line-of-sight ( Muxlow fc Garrington 199 1|) . This inclination readily 



explains the seemingly large rotation of the jet from the axis defined by N and S to that 
defined by C and E. 

The performance of the VLA at A18cm has improved markedly since 1980, and new 
observations should detect or exclude these features with high significance. We are aware 
of a very deep VLA observation ( [Harvanek et at. 1996j ) at A18cm and A3. 6cm; the longer 



wavelength data should be able to confirm GE, GN, and GNE, and if GE is detected at 
A3. 6cm it may be possible to determine whether the mass model is singular, or whether GE 
consists of two merging images. 
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09 57 59 58 57 56 

RIGHT ASCENSION 



Fig. 1. — Contour plot of A18cm A array map of 0957+561 on 1980 December 16. The 
cross hair marks the position at which the quasar B component has been subtracted from 
the map. The source just to the north of B is G, the lensing galaxy. Contour levels are 
-0.10%, 0.10%, 0.20%, 0.28%, 0.40%, 0.57%, 0.80%, 1.13%, 1.60%, 2.26%, 3.2%, 4.53%, 
6.40%, 9.05%, 12.8%, 18.1%, 25.6%, 36.2%, and 51.2% of the peak intensity of 181 mJy 
beam -1 . The noise level is 105^Jy beam -1 . The box in the lower left shows the beam 
FWHM ellipse. 
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Fig. 2. — Contour plot of A6cm map of 0957+561 from co-added observations in A, B, and 
C arrays (data sets #4 - #8, Table 1). The cross hairs (+) are the positions from which 
models of the components B (to the south) and G were subtracted. The crosses ( x ) are, east 
to west, the positions of HST components "blob 2" and "blob 3." The circles are positions 
along the HST arc, the outer two being the approximate extent and the inner two being 
"knot 1" and "knot 2." Contour levels are -0.25%, 0.25%, 0.35%, 0.50%, 0.63%, 0.75%, 
0.88%, 1.00%, 1.13%, 1.60%, 2.26%, 3.20%, 4.53%, 6.40%, 9.05%, 12.8%, and 51.2% of the 
peak intensity of 41.8 mJy beam -1 . The noise level is 39/xJy beam" 1 . The box in the lower 
left shows the beam FWHM ellipse. 
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Fig. 3. — Lens model constrained to the HST components, showing disposition of radio 
components. The HST and radio components are shown as dark and light contours, 
respectively. The source plane shows how the source would appear without lensing; the 
caustics separate regions of multiple imaging. The image plane shows the model source seen 
through our lens model; the critical lines divide the images. The location of G is shown on 
the image plane, at the center of the lens model. Note that the HST arc is probably formed 
by the eastern end of the source that yields the HST blobs. 
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Table 1. Archival VLA Data Presented in this Letter. 



N° 


Obs. Date 


A (cm) 


Duration (hrs) a 


VLA Coring. 


1 


16DEC80 


18 


3.0 


A 


2 


21AUG90 


6 


0.8 


B 


3 


17MAR88 


6 


0.9 


C 


4 


27SEP87 


6 


1.5 


A 


5 


20JUL87 


6 


0.8 


A 


6 


09DEC87 


6 


1.1 


B 


7 


02JUN85 


6 


1.2 


B 


8 


11JAN87 


6 


1.6 


C 



a The duration is the effective VLA observing time devoted to 
0957+561, excluding time spent on calibration and instrumental 
difficulties. 
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Table 2. Faint Emission Features Described in this Letter 



Feature 




// n 
(Aa ,A<5 ) a 


A6cm Peak b 


A6cm Flux 


A18cm Flux c 








mJy beam -1 


mJy 


mJy 


N 


(-4.5 , 


18.8 ) ± (0.5 ,0.5 ) 


0.440±0.053 d 


0.840±0.110 d 


10.15±0.66 


S 


( 1-5 ,- 


-11.7 ) ± (0.5 ,0.5 ) 


0.270±0.053 d 


1.030±0.149 d 


8.87±0.61 


Rl 


( 5.38, 


0.80) ± (0.59,1.17) 


0.366±0.043 c 


3.646±0.181 c 


6.61±0.31 


R2 


(-1.28, 


-0.80) ± (0.26,0.21) 


0.250±0.039 f 


2.270±0.159 f 




GE 


( 1-24, 


0.68) ± (0.04,0.04) 


0.486±0.039 f 


2.130±0.093 f 




GN 


(-0.24, 


1.68) ± (0.12,0.24) 


0.373±0.039 f 


0.930±0.066 f 




GNE 


( 2.24, 


2.08) ± (0.34,0.22) 


0.357±0.039 f 


1.186±0.073 f 





Positions relative to B, a=09 h 57 m 57 s A2±(fm, £=56°08'16"40±0"1 (B1950) 

b Flux uncertainties are based on the measured map noise away from source emission. 
For the components GE, GN, and GNE, the error is dominated by the deconvolution 
algorithm, and the quoted errors are likely underestimated. 

c A18cm, beam FWHM ellipse 1.94" x 1.47", map rms 105/zJy beam -1 
d A6cm, beam FWHM ellipse 3.51" x 2.57", map rms 53/iJy beam -1 
c A6cm, beam FWHM ellipse 1.56"xl.37", map rms 43;uJy beam -1 
f A6cm, beam FWHM ellipse 0.75"x0.69", map rms 39^Jy beam -1 
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Table 3. Best Fit Lens Model Parameters 



Parameter Best Fit Value 



Aa a +0"181±0"001 

A<5 a +1'.'019±0'.'001 

b 2'.'88±0'.'021 

9 C (<0'/02) 

P 0.88±°;°2 

e 0.538±0.018 

(j) h +65.2±0.5 

Q 2 1.03±0.02 



a offsets in right ascension and declination from the B quasar image. 
b orientation in degrees from north through east 



